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Abstract 



The D0 collaboration reports on a search for the Standard Model top quark 
in pp collisions at yfs = 1.8 TeV at the Fermilab Tevatron, with an integrated 
luminosity of approximately 50 pb _1 . We have searched for ti production in 
the dilepton and single- lepton decay channels, with and without tagging of 
6-quark jets. We observed 17 events with an expected background of 3.8 ± 0.6 
events. The probability for an upward fluctuation of the background to pro- 
duce the observed signal is 2 x 10~ 6 (equivalent to 4.6 standard deviations). 
The kinematic properties of the excess events are consistent with top quark 
decay. We conclude that we have observed the top quark and measure its 
mass to be 19913,1 (stat.) ±22 (syst.) GeV/c 2 and its production cross section 
to be 6.4 ±2.2 pb. 
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In the Standard Model (SM), the top quark is the weak isospin partner of the b quark. 
The D0 collaboration published a lower limit on the mass of the top quark of 131 GeV/c 2 , 
at a confidence level (CL) of 95%, based on an integrated luminosity of 13.5 pb -1 M. A 
subsequent publication reported the top quark production cross section as a function of 
the assumed top quark mass. In that analysis, we found nine events with an expected back- 
ground of 3.8 ± 0.9 events (statistical significance 1.9 standard deviations) corresponding to 
a production cross section of 8.2 ± 5.1 pb for an assumed top quark mass of 180 GeV/c 2 . 
The CDF collaboration published evidence for top quark production with a statistical signif- 
icance of 2.8 standard deviations, a top quark of mass 174±10t]J GeV/c 2 , and a production 
cross section of 13.9jj|;g pb ||. Precision electroweak measurements predict a SM top quark 
mass of approximately 150-210 GeV/c 2 , depending on the mass of the Higgs boson ||. In 
the present paper, we report new results from the D0 experiment that firmly establish the 
existence of the top quark. 

We assume that the top quark is pair-produced and decays according to the minimal 
SM (i.e. ti — > W + W~hh). We have searched for the top quark in channels where both W 
bosons decayed leptonically (e// + jets, ee + jets, and //// + jets) and in channels where just 
one W boson decayed leptonically (e + jets and // + jets). The single-lepton channels were 
subdivided into 6-tagged and untagged channels according to whether or not a muon was 
observed consistent with b — > // + Jf. The muon-tagged channels are denoted e+jets/// and 
// + jets///. 

Here we present an analysis based on data collected at the Fermilab Tevatron at ^/s = 
1.8 TeV with an integrated luminosity of 44-56 pb -1 , depending on the channel. In the 
present analysis, the signal-to-background ratio for a high mass top quark was substantially 
improved relative to Ref. 0. An optimization of the selection criteria was carried out using 
Monte Carlo top quark events for signal and our standard background estimates. The result 
of this procedure was a factor of 3.7 better background rejection while retaining 70% of the 
acceptance for 180 GeV/c 2 top quarks. This corresponds to a signal-to-background ratio of 
1:1 for a top quark mass of 200 GeV/c 2 , assuming the expected SM top cross section ||. 
The improved rejection arises primarily by requiring events to have a larger total transverse 
energy. 

The D0 detector and data collection systems are described in Ref. ||. The triggers and 
reconstruction algorithms for jets, electrons, muons, and neutrinos were the same as those 
used in our previous top quark searches |l|Q]. 

The signature for the dilepton channels was defined as two isolated leptons, at least two 
jets, and large missing transverse energy Ifa. The signature for the single-lepton channels 
was defined as one isolated lepton, large I$t, and a minimum of three jets (with muon tag) 
or four jets (without tag). The minimum transverse momentum px of tagging muons was 
4 GeV/c. Requirements pertaining to the magnitude and direction of the I^t, the aplanarity 
of the jets A, and the allowed ranges of pseudorapidity r\ were similar to Ref. 0. Muons 
were restricted to \rj\ < 1 for the last 70% of the data because of forward muon chamber 
aging. Events in the //// + jets and // + jets/// channels were required to be inconsistent 
with the Z + jets hypothesis, based on a global kinematic fit. The principal difference 
between the present analysis and the analysis of Ref. was the imposition of a minimum 
requirement in all channels on a quantity H?, which we defined as the scalar sum of the 
transverse energies E T of the jets (for the single-lepton and //// + jets channels) or the 
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scalar sum of the E^s of the leading electron and the jets (for the e/i + jets and ee + jets 
channels). The kinematic requirements for our standard event selection for all seven channels 
are summarized in Table |I[ In addition to the standard selection, we defined a set of loose 
event selection requirements, which differed from the standard set by the removal of the Ht 
requirement and by the relaxation of the aplanarity requirement for e + jets and /i + jets 
from A > 0.05 to A > 0.03. 

For the dilepton channels, the main backgrounds were from Z and continuum Drell-Yan 
production (Z, 7* — > ee,/i/i, and tt), vector boson pairs (WW 7 , WZ), heavy flavor (bb and 
cc) production, and backgrounds with jets misidentified as leptons. For the single-lepton 
channels, the main backgrounds were from jy + jets, Z + jets, and multijet production with 
a jet misidentified as a lepton. The method for estimating these backgrounds was the same 
as in our previously published analyses 

Ht is a powerful discriminator between background and high-mass top quark production. 
Figure |T] shows a comparison of the shapes of the Ht distributions expected from background 
and 200 GeV/c 2 top quarks in the channels (a) e/i + jets and (b) untagged single-lepton + 
jets. We have tested our understanding of background Ht distributions by comparing data 
and calculated background in background-dominated channels such as electron + two jets 
and electron + three jets (Fig. |2|). The observed Ht distribution agrees with the background 
calculation, which includes contributions from both W + jets as calculated by the vecbos 
Monte Carlo [0] and multijet events. 

The acceptance for tt events was calculated using the isajet event generator || and 
a detector simulation based on the geant program ||. As a check, the acceptance was 
also calculated using the herwig event generator |ID |. The difference between isajet and 
herwig was included in the systematic error. 

From all seven channels, we observed 17 events with an expected background of 3.8 ±0.6 
events (see Table [D]). Our measured cross section as a function of the top quark mass 
hypothesis is shown in Fig. |3|. Assuming a top quark mass of 200 GeV/c 2 , the production 
cross section is 6.3±2.2 pb. The error in the cross section includes an overall 12% uncertainty 
in the luminosity. The probability of an upward fluctuation of the background to 17 or more 
events is 2 x 10~ 6 , which corresponds to 4.6 standard deviations for a Gaussian probability 
distribution. We have calculated the probability for our observed distribution of excess 
events among the seven channels and find that our results are consistent with top quark 
branching fractions at the 53% CL. Thus, we observe a statistically significant excess of 
events and the distribution of events among the seven channels is consistent with top quark 
production. 

Additional confirmation that our observed excess contains a high-mass object comes 
from the invariant masses of jet combinations in single-lepton + jets events. For this anal- 
ysis, we selected single-lepton + four-jet events using the loose event selection require- 
ments (27 events). An invariant mass analysis was performed, based on the hypothesis 
tt — > W + W~bb — > ivqqbb. One jet was assigned to the semileptonically decaying top quark 
and three jets were assigned to the hadronically decaying top quark. The jet assignment al- 
gorithm attempted to assign one of the two highest Et jets to the semileptonically decaying 
top quark and to minimize the difference between the masses of the two top quarks. The 
invariant mass of the three jets assigned to the hadronically decaying top quark is denoted 
by m 3 j. The invariant mass of the pair of hadronically decaying top quark jets with the 
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smallest invariant mass is denoted by m^. Figure ^ shows the distribution of vs. m2j for 
(a) background (W +jets and multijet) (b) 200 GeV/c 2 top Monte Carlo, and (c) data. The 
data are peaked at higher invariant mass, in both dimensions, than the background. Based 
only on the shapes of the distributions, the hypothesis that the data are a combination of 
top quark and background events (60% CL) is favored over the pure background hypothesis 
(3%CL). 

To measure the top quark mass, single-lepton + four-jet events were subjected to 2- 
constraint kinematic fits to the hypothesis ti — > W + W~bb — > iuqqbb. Kinematic fits were 
performed on all permutations of the jet assignments of the four highest jets, with the 
provision that muon-tagged jets were always assigned to a 6-quark in the fit. A maximum 
of three permutations with x 2 < 7 (two degrees of freedom) were retained, and a single 
X 2 -probability- weighted average mass ("fitted mass") was calculated for each event. Monte 
Carlo studies using the isajet and herwig event generators showed that the fitted mass 
was strongly correlated with the top quark mass. Gluon radiation, jet assignment combina- 
torics, and the event selection procedure introduced a shift in the fitted mass (approximately 
—20 GeV/c 2 for 200 GeV/c 2 top quarks), which was taken into account in the final mass 
determination. 

Eleven of the 14 single-lepton + jets candidate events selected using the standard cuts 
were fitted successfully. Figure |5](a) shows the fitted mass distribution. An unbinned like- 
lihood fit, incorporating top quark and background contributions, with the top quark mass 
allowed to vary, was performed on the fitted mass distribution. The top quark contribution 
was modeled using isajet. The background contributions were constrained to be consis- 
tent with our background estimates. The likelihood fit yielded a top quark mass of 199^25 
(stat.) GeV/c 2 and described the data well. 

To increase the statistics available for the mass fit, and to remove any bias from the 
standard Ht requirement, we repeated the mass analysis on events selected using the loose 
requirements. Of 27 single-lepton + four-jet events, 24 were fitted successfully. The removal 
of the Ht requirement introduced a substantial background contribution at lower mass in 
addition to the top signal, as shown in Fig. 5(b). A likelihood fit to the mass distribution 
resulted in a top quark mass of 199^2? (stat.) GeV/c 2 , consistent with the result obtained 
from the standard event selection. The result of the likelihood fit did not depend significantly 
on whether the normalization of the background was constrained. Using herwig to model 
the top quark contribution resulted in a mass 4 GeV/c 2 below that found using isajet. 
This effect was included in the systematic error. The total systematic error in the top quark 
mass is 22 GeV/c 2 , which is dominated by the uncertainty in the jet energy scale. 

In conclusion, we report the observation of the top quark. We measure the top quark 
mass to be 1991 J? (stat.) ±22 (syst.) GeV/c 2 and measure a production cross section of 
6.4 ± 2.2 pb at our central mass. 

We thank the Fermilab Accelerator, Computing, and Research Divisions, and the support 
staffs at the collaborating institutions for their contributions to the success of this work. 
We also acknowledge the support of the U.S. Department of Energy, the U.S. National 
Science Foundation, the Commissariat a L'Energie Atomique in France, the Ministry for 
Atomic Energy and the Ministry of Science and Technology Policy in Russia, CNPq in 
Brazil, the Departments of Atomic Energy and Science and Education in India, Colciencias 
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FIGURES 




H T (GeV) 



FIG. 1. Shape of Ht distributions expected for the principal backgrounds (dashed line) and 
200 GeV/c 2 top quarks (solid line) for (a) e/x + jets and (b) untagged single-lepton + jets. 




FIG. 2. Observed distributions (points) compared to the distributions expected from back- 
ground (line) for rfj< > 25 GeV/c and (a) e + > 2 jets and (b) e + > 3 jets. 
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FIG. 3. D0 measured ti production cross section (solid line with one standard deviation error 
band) as a function of assumed top quark mass. Also shown is the theoretical cross section curve 
(dashed line) ||. 
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(a) (b) (c) 




FIG. 4. Single-lepton + jets two-jet vs. three-jet invariant mass distribution for (a) background, 
(b) 200 GeV/c 2 top Monte Carlo (isajet), and (c) data. 




100 200 100 200 

Fitted Mass (GeV/c 2 ) 

FIG. 5. Fitted mass distribution for candidate events (histogram) with the expected mass 
distribution for 199 GeV/c 2 top quark events (dotted curve), background (dashed curve), and the 
sum of top and background (solid curve) for (a) standard and (b) loose event selection. 
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TABLES 



TABLE I. Minimum kinematic requirements for the standard event selection (energy in GeV). 





Leptons 


Jets 








Channel 


£r(e) 


Prfji) 




Ex 


■f i 




A 


e// + jets 


15 


12 


2 


15 


20 


120 




ee + jets 


20 




2 


15 


25 


120 




//// + jets 




15 


2 


15 




100 




e + jets 


20 




4 


15 


25 


200 


0.05 


fj, + jets 




15 


4 


15 


20 


200 


0.05 


e + jets/// 


20 




3 


20 


20 


140 




\x + jets/// 




15 


3 


20 


20 


140 





TABLE II. Efficiency x branching fraction (e x B) using standard event selection and the 
expected number of top quark events ({N)) in the seven channels, based on the central theoretical 
ti production cross section of Ref. ||, for four top masses. Also given are the expected background, 
integrated luminosity, and the number of observed events in each channel. 



m t (GeV/c 2 ) 


e/i + jets 


ee + jets 


fifi + jets 


e + jets 


fi + jets 


e + jets/ fi 


fi + jets//i 


ALL 


e x B(%) 


0.17 ±0.02 


0.11 ±0.02 


0.06 ±0.01 


0.50 ±0.10 


0.33 ±0.08 


0.36 ±0.07 


0.20 ±0.05 




140 (N) 


1.36 ± 0.21 


1.04 ±0.19 


0.46 ±0.08 


4.05 ±0.94 


2.47 ±0.68 


2.93 ±0.68 


1.48 ± 0.42 


13.80 ±2.07 


e x B{%) 


0.24 ±0.02 


0.15 ±0.02 


0.09 ±0.02 


0.80 ±0.10 


0.57 ±0.13 


0.50 ±0.08 


0.25 ± 0.06 




160 (N) 


0.94 ±0.13 


0.69 ±0.12 


0.34 ±0.07 


3.13 ±0.54 


2.04 ±0.53 


1.95 ±0.39 


0.92 ± 0.24 


10.01 ± 1.41 


e x S(%) 


0.28 ± 0.02 


0.17 ±0.02 


0.10 ±0.02 


1.20 ±0.30 


0.76 ±0.17 


0.56 ±0.09 


0.35 ± 0.08 




180 (N) 


0.57 ±0.07 


0.40 ± 0.07 


0.19 ±0.04 


2.42 ±0.67 


1.41 ±0.36 


1.14 ±0.22 


0.64 ±0.16 


6.77 ± 1.09 


e x B{%) 


0.31 ± 0.02 


0.20 ±0.03 


0.11 ±0.02 


1.70 ±0.20 


0.96 ±0.21 


0.74 ±0.11 


0.41 ± 0.08 




200 (N) 


0.34 ± 0.04 


0.25 ±0.05 


0.11 ±0.02 


1.84 ±0.31 


0.95 ±0.24 


0.81 ±0.16 


0.41 ± 0.10 


4.71 ±0.66 


Background 


0.12 ± 0.03 


0.28 ±0.14 


0.25 ±0.04 


1.22 ±0.42 


0.71 ±0.28 


0.85 ±0.14 


0.36 ± 0.08 


3.79 ±0.55 


/' Cdt (pb- 1 ) 


47.9 ±5.7 


55.7 ±6.7 


44.2 ±5.3 


47.9 ±5.7 


44.2 ±5.3 


47.9 ±5.7 


44.2 ± 5.3 




Data 


2 





1 


5 


3 


3 


3 


17 



12 



